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1.  Introduction 


The  U.S.  Army  Research  Laboratory  is  embarking  on  a  program  of  holographic  interferometry  to 
understand  its  limitations  for  the  purpose  of  remote  sensing.  A  pre-requisite  for  this  work  is  to 
record  a  digital  hologram.  When  two  holograms  are  recorded  and  processed,  using  two  different 
wavelengths,  three-dimensional  images  are  formed.  Such  images,  depending  on  the  wavelengths 
selected,  show  depth  information  that  can  approach  microscopic  dimensions.  Other  advantages 
of  using  holographic  techniques  are  the  large  depth  of  field,  the  lack  of  a  need  for  mechanical 
focusing  mechanisms,  and  the  perfect  image  reconstruction  that  possess  both  phase  and 
amplitude  information  of  the  object  being  examined,  instead  of  just  the  intensity  information  that 
is  in  a  regular  photograph  ( 1 ). 

Although  this  work  can  and  has  been  performed  with  chemical  films,  digital  capture  and 
processing  of  images  can  be  accomplished  much  more  quickly.  For  this  reason,  mathematical 
methods  that  can  be  applied  digitally  have  been  developed  to  process  the  hologram  interference 
patterns  into  images.  Of  particular  interest  are  methods  that  use  Fourier  transforms,  because  this 
type  of  transform  is  well  known  and  has  been  implemented  in  many  image  processing  programs, 
such  as  IDL  and  MATLAB,  just  to  name  a  few. 


2.  Methods,  Assumptions,  and  Procedure 


2.1  Mathematical  Analysis,  Part  1 — Establishing  a  Mathematical  Framework: 

Developing  the  Helmholtz-Kirchhoff  Equation 

The  purpose  of  the  following  analysis  is  to  develop  a  method  to  transform  a  hologram,  which  is 
an  interference  pattern  of  light  obtained  in  the  laboratory,  into  an  image.  There  are  two  main 
ways  to  do  this.  The  first  method,  which  is  described  below,  uses  a  Fourier  transform.  The 
second  method,  not  discussed  here,  uses  convolution.  We  chose  to  use  the  Fourier  transform, 
because  it  has  been  optimized  in  many  types  of  computer  applications  that  we  use,  and  because 
signal  processing  hardware  exists  that  can  optimize  the  process  of  performing  the  transform.  We 
show  and  discuss  this  analysis  from  a  pedagogical  point  of  view,  making  the  connecting  steps  in 
the  analysis  clear  and  discussing  what  the  final  results  of  the  analysis  mean. 

The  analysis  begins  with  Stake’s  theorem  equation  1,  which  relates  a  volume  integral  to  a 
surface  integral.  Figure  1  shows  the  geometry  associated  with  equation  1. 
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Figure  1 .  The  surface,  volume,  and  normal  vector 

that  relate  the  two  integrals  in  Stoke ’s  theorem. 


HI (UlA2U2  - U2A2Ul )dv  =  || - U2AUl )  •  ft  ds  (l) 

V  5 

In  this  case,  the  functions  Ui  and  U2  in  equation  1  are  light  waves.  This  means  that  these 
functions  are  solutions  to  the  wave  equation.  Maxwell’s  equations  show  that  the  volume  integral 
on  the  left  hand  side  of  equation  1  goes  to  zero.  Equation  1  then  reduces  to 

0  =  || (Ux AU2  -  U2AUl )»nds  (2) 

5 

An  analogy  with  the  scaler  dot  product  ( A  •  B )  is  useful  in  further  evaluating  equation  2.  The 

vector  product  of  A  and  B  means  that  the  component  of  vector  A  is  projected  in  the  B  direction. 
(A U  •  ft)  is  therefore  the  change  in  Ui  in  the  n  direction,  so  equation  2  can  be  written  as 

0=l\(u>^-u^)ds  O) 

on  on 


Consider  now  figure  2  to  proceed  in  the  analysis. 


Figure  2.  Two  surfaces,  one  of  which  surrounds 
a  point,  used  to  evaluate  equation  3. 
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The  complex  amplitude  on  Si  is  used  to  calculate  the  complex  amplitude  on  S2  surrounding  the 
small  sphere,  P.  The  unknown  wavefront  at  S2,  which  is  a  complex  amplitude,  is  evaluated  as  a 
spherical  wavefront.  This  analysis  works  well  for  holography,  where  the  individual  point  “P” 
might  be  a  pixel  in  a  charge-coupled  device  (CCD)  array  or  the  emulsion  of  a  film.  The 
unknown  wave  at  S2  has  the  form 


(4) 


where  ‘W’  is  the  distance  from  point  P,  and  “k”  is  the  wave  number.  The  integral  is  now 
evaluated  over  each  surface  using  equation  4,  where  r2  becomes  infinitesimally  small. 
Substitution  of  equation  4  into  equation  3  over  Si  becomes 
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The  surface  integral  at  surface  S2  is  evaluated  using  the  fact  that  the  radius  r2  is  in  the  same 
direction  as  no. 

d  _  d 

dn2  dr2 


After  differentiating,  the  surface  integral  for  surface  S2  becomes 
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Since  r2  (from  figure  2)  is  reduced  to  zero,  r2=s,  and  dS2=s"dQ,  equation  reduces  to 
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which,  when  integrated,  yields 
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Terms  1  and  3  in  equation  9  go  to  zero  as  s  goes  to  zero.  Term  2  and  equation  9  then  reduce  to 
4tiUi  as  s  goes  to  zero. 
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Since  the  complex  amplitudes  on  the  two  surfaces  (Si  and  S2  in  figure  2)  are  equal,  equation  5 
and  equation  7  are  equal.  Equation  7  reduces  to  4jiUi,  as  discussed  previously,  so  the  equality 
between  Si  and  S2  becomes 
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(10) 


For  the  digital  holography  application  being  considered,  think  of  Up  as  being  the  amplitude  at  the 
pixel  of  a  CCD.  This  result  is  called  the  Helmholtz-Kirchhoff  equation.  Equation  10  is  difficult 
to  solve  analytically  and  time  consuming  if  approached  computationally. 

2.2  Mathematical  Analysis,  Part  2 — Integrating  an  Aperture  Instead  of  a  Full  Sphere: 

The  Kirchhoff  Assumptions 

Surface  Si  of  figure  2  is  a  closed  surface  encompassing  a  full  hemisphere.  The  Kirchhoff 
assumption  is  that  the  full  sphere  is  not  illuminated;  only  a  small  portion  (an  aperture)  of  that 
sphere  is  illuminated.  The  following  assumptions  then  apply: 

U]  =  0  and  ^^-  =  0  (11) 

on 


on  the  surface  Si  except  at  the  aperture.  The  following  geometry  is  useful  for  visualizing  the 
situation  and  understanding  trigonometric  relations  used  in  the  approximation. 

The  aperture,  which  replaces  Si,  is  approximated  as  planar  because  it  is  such  a  small  portion  of 
the  surface  Si.  Using  the  coordinates  in  figure  3,  the  spherical  wavefront  of  the  source  is 

„  Aelkr' 
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Figure  3.  Geometries  of  the  aperture  that  replaces  Si. 


The  normal  to  the  aperture  does  not  change  (as  it  does  on  the  surface  Si)  and  is  the  z-axis,  so 
equation  12  can  be  differentiated  as  follows: 
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Light  is  coming  to  the  aperture  from  the  source  of  light  and  from  the  observing  point,  which  is 
the  object  that  is  illuminated  to  create  a  hologram.  Light  from  both  of  these  sources  reflect  back 
and  interfere  at  the  aperture.  The  light  from  the  observing  point  is  LL  and 
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A  final  assumption  is  made,  and  that  is  that  the  source  and  object  are  more  than  1  cm  from  the 
aperture,  since  k=27i/k~  100,000  for  visible  light,  and  in  equations  13  and  14,  if  ri  and  r2  are  on 
the  order  of  1  cm  or  greater,  then  1/ri  and  l/r2  can  be  neglected  with  respect  to  k.  This  means 
that 


l 

rl 


(15) 


Equation  10,  which  is  copied  below  as  equation  15,  is  the  starting  point  for  applying  the  aperture 
assumptions. 


The  numbers  in  parenthesis  are  the  equations  are  used  to  get  from  equation  16  to 
Combining  like  terms  reduced  equation  17  to 

U P  =  — — —  IT  -T—  etk ( 'i  +,2 )  ( cos  ^  -  cos  <p2 )  dxdy 

Aperture  1  L 

Since  k=2jt/A,,  equation  1 8  can  be  further  reduced  to 

UP  =  — JJ  _L  elk(n+r2)  ^CQS  ^  _  CQS  ^  ^  fady 

Aperture  ^1^2 


equation  17. 

(18) 

(19) 


This  result  is  called  the  Kirchhoff  diffraction  integral  and  is  still  difficult  to  integrate  in  most 
situations.  The  reason  is  that  as  the  integration  element,  ds=dxdy,  moves  around  the  aperture  in 
figure  3  and  the  distances  ri  and  r2  along  with  the  direction  cosines  vary. 

2.3  Mathematical  Analysis,  Part  3 — The  Fresnel  Approximation 

Figure  3,  which  is  a  global  coordinate  system  in  x-y-z  variables,  is  modified  by  adding  a  local 
coordinate  system  in  the  single  dimension  of  cj.  This  is  shown  in  figure  4. 
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The  Fresnel  approximation  further  assumes  that  points  Q  and  P  are  located  a  distance  from  the 
aperture  that  is  larger  than  the  aperture  size.  This  means  that  the  points  are  “far”  away  and  the 
aperture  is  “small.”  The  factor  l/r^  that  appears  in  the  integral  of  equation  19  can  now  be  taken 
outside  of  the  integral  and  replace  by  l/riW.  The  primed  values  show  that  the  values  are 
measured  with  respect  to  the  local  coordinate  system.  Since  the  points  are  “far”  from  the 
aperture,  the  directions  cosines  not  vary  much  and  can  be  considered  a  constant,  so  the 
trigonometric  equality  may  be  used  to  replace  (cosqq-coscpi)  with  2coscp.  These  approximations 
remove  certain  terms  from  equation  19,  because  they  are  constant.  However,  the  values  in  the 
exponent  of  equation  19  cannot  be  treated  as  constants,  because  the  exponent  oscillates  as  the 
integration  element  moves  over  the  aperture.  With  these  assumptions  and  considerations  in 
mind,  equation  19  is  simplified  to 

TT  COS  (j)  |*  ik(r1+r2)  j  z 

=  — —  J  <>  -W  (20) 

1  2  Aperture 

The  Pythagorean  and  binomial  theorems  are  now  used  to  simplify  the  integrand  (2). 
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(21) 


Because  points  Q  and  P  are  far  from  the  aperture,  the  three  terms  in  the  parenthesis  of 
equation  21  are  small,  and  the  binomial  theorem  may  be  used.  Equation  21  then  becomes 

xf  £2  X& 

r\  =  z,  +  —  +  - - — 

2zj  2z,  z. 

Substituting  equation  22  and  the  corresponding  value  for  n  into  equation  20  yield 


ik(  x?  x? 
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Aperture 


This  result  is  called  the  Fresnel  approximation  to  the  Kirchhoff  diffraction  integral.  A  final 
simplification  is  required  to  make  the  result  useful  for  most  practical  applications. 

2.4  Mathematical  Analysis,  Part  4 — The  Fraunhofer  Approximation 

The  central  point  of  the  Fraunhofer  approximation  is  to  make  the  first  exponential  term  in 
equation  23  nearly  equal  to  one.  This  is  the  case  if  t,  is  small  compared  to  the  distances  of  the 
Illumination  source,  Q,  and  the  illuminated  target,  P,  which  are  at  distances  of  zi  and  Z2, 
respectively.  This  assumptions  change  equation  23  to 

‘±U+4)  -ik{*AM  1 

2  l.  Z2  )  f  e  U'l  Z2  ) 

Aperture 

Y - ^ - Y - '  (24) 

Illumination  Spherical  Fourier 

at  Wave  Transform 

Aperture 


U  _ -At  cos <f>  jk(Zl+z9, 
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3.  Conclusion 


Equation  24  is  the  basis  upon  which  ARL’s  holographic  interferometry  efforts  are  built,  and  this 
is  the  reason  such  a  long  discourse  is  devoted  to  the  subject.  The  final  result  of  equation  24 
shows  that  diffraction  at  an  aperture  is  really  a  Fourier  transforming  process.  The  optical 
information  of  the  object  at  point  P  has  been  transformed  into  the  frequency  information  that 
comprises  a  hologram,  in  holography,  at  the  aperture.  The  aperture  may  be  a  film  or  a  digital 
recording  device. 
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